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The Hamiltonian Mean Field (HMF) model is a paradig- 
matic toy model for long-range Hamiltonian system, which 
has raised recently a great attention and has been investi- 
gated in detail, see for example [2] and references therein. 
We have recently published two papers on this model [3,4] 
concerning the possible application of a generalized version 
of the Central Limit Theorem (CLT) in accordance with 
g-statistics [5-8]. In these papers, following ref. [7], we 
constructed probability density functions (Pdf) of quan- 
tities ?/'s expressed as a finite sum of n stochastic vari- 
ables, selected along the deterministics time evolutions of 
tjhe N rotors at fixed time intervals 6. We found that, 
along the metastable quasistationary states (QSS) of the 
HMF model , g-Gaussian-like Pdfs appear to emerge as 
attractors [3,4]. In ref. [1] the authors advance several 
misleading claims suggesting that the results found in pa- 
pers [3,4] are erroneus and that no q-Gaussian attractors 
exist for correlated variables extracted from HMF dynam- 
ics. Two crucial points raised by Figuereido et al. are 
the following: (i) by doing the calculations with more ac- 
curacy and better statistics, it would be possible to show 
that a q-Gaussian attractor does not exist, since there is 
an evident cut-off in the tails of the final distributions; 
(ii) by considering increasingly large size systems, the Pdf 
would gradually converge to a Gaussian. The point (i), 
which is illustrated in fig.l of [1], can be easily rebutted 
by considering that, like the usual CLT, the g-CLT strictly 
applies only if n goes to infinity. In the example shown 
in [1] only n = 50 (with N = 100 and S = 40) was con- 
sidered, therefore one expects deviations in the tails, and 
elsewhere. Here we present new calculations for the case 
N = 20000, which indicate that this phenomenon is true 
whatever is the shape of the final attractor. It is also im- 
portant to remind that, in ref. [4] (ref. [33] of [1]), we 



observed the existence of three classes of trajectories for 
the same water-bag (Ml) initial conditions, each one of 
them showing a different central limit behavior: class 1, 
with a g-Gaussian-like attractor, class 3 with a Gaussian 
attractor and class 2 with an intermediate behavior. In the 
upper row of fig.l we show the Pdf of a typical event of 
class 1 for three increasing values of n (n = 200, n = 1000 
and n = 2000), while in the middle row the correspond- 
ing behavior is illustrated for a typical event of class 3, 
for the same values of n. Being in all the simulations 
S = 100, the sums have been performed, respectively, for 
20000, 100000 and 200000 time steps. The time evolution 
for the temperature in both cases is plotted in the bottom 
panel. Notice that the cut-off in the tails of the distri- 
butions, clearly visible for n = 200, disappears slowly by 
increasing the number of summands n, in both cases. The 
convergence starts, for small n, from the central part of 
the Pdf and then extends to the tails, producing in one 
case (class 3) a Gaussian-like and, in the other case (class 
1), a g-Gaussian-like curve as a final attractor. The other 
main point addressed by fig.l is that at least two kinds of 
attractors coexist in the QSS regime and their occurrence, 
as discussed in ref. [4] , depends very sensitively on the dif- 
ferent realizations of the initial conditions. The authors 
of ref. [1], although quoting our paper, fail in taking into 
account this fundamental point: in their fig. 3 they con- 
sider a system of N = 100000 rotators and plot a Pdf 
for 50 events without any reference to their corresponding 
class, and it is evident that the trajectories used belonged 
to a random mixture of classes 1, 2 and 3 (see [4]). Ac- 
tually, mixing events or not distinguishing between events 
of different classes leads to unreliable conclusions, since 
the dynamics and the correlations of the various trajec- 
tories can be very different. We notice that although the 
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Fig. 1: Time evolution of two events with the same Ml initial 
conditions, more precisely with N = 20000 at U = 0.69, but 
belonging to different classes (1 and 3). The evolution towards 
the final attractor appears to be a Gaussian-like for the event 
of class 3 and a q-Gaussian-like for the event of class 1. The q- 
Gaussian is a generalization of Gaussian which emerges in the 
context of nonextensive statistical mechanics [5] and is defined 
asG„(x) = A(l-(l~q)f3x 2 ) 1/1 - q . Here we have q = 1.42±0.1, 
/3 = 1.3 ± 0.1 and A = 0.55, and the resulting curve is plotted 
in all the panels toghether with a standard Gaussian. See text 
for further details. 



events of class 1 seem to diminish with N, those of class 
2 increase more than those of class 3 [4] . Finally, in order 
to address point (ii) of the criticism raised in ref. [1], we 
illustrate in fig. 2 how the attractor shown in the upper 
row of fig. 1 evolves its shape by increasing both the size of 
the system N and, consistently, the number of summands 
n, but leaving constant the ratio N/n. This latter choice 
is however not a necessary condition, here it is done only 
because the QSS lifetime increases with N. The calcu- 
lations refer here to three typical events of class 1 with 
sizes N = 10000, 50000 and 100000 (note that the range 
of the ordinate was enlarged with respect to fig. 1). The 
gradual filling of the tails by increasing N strongly sug- 
gests that, for N and n going to infinity, a g-Gaussian-like 
asymptotic shape might well be eventually reached, thus 
reinforcing its role of real attractor for the QSS regime. 
In conclusion, the evidence presented in [1] by no means 
invalidates the occurrence of the g-generalizcd CLT for 
the HMF system and in general for long-range Hamilto- 
nian systems. On the other hand, if these systems are 
non-ergodic (as the authors themselves admit and as we 
believe we have clearly established [3,4]), it also remains 
unclear why should one not expect, for the time-averaged 
distribution, a non-Gaussian attractor in the QSS regime. 
In summary, we presented further numerical evidence that 
standard statistical mechanics has severe problems in ex- 
plaining the anomalous long-standing QSS regime of the 




Fig. 2: For U = 0.69 we show the Pdfs obtained considering 
single events of class 1 for different sizes N, with S = 100 
and n = N/10. Again the indications for a g-Gaussian-like 
attractor becomes stronger and stronger when sending both N 
and n to infinity. See text for further details. 



HMF model, showing that in this respect g-statistics re- 
mains a very good candidate. 

REFERENCES 

[1] FlGUEIREDO A., ROCHA FlLHO T.M., AMATO M.A., Eu- 

roph. Lett, 83 (2008) 30011 
[2] Dauxois T., Latora V. , Rapisarda A. , Ruffo S. and 
Torcini A., Lecture Notes in Physics, edited by T. Daux- 
ois, S. Ruffo, E. Arimondo, M. Wilkens, Vol. 602 
2002, p. 458; Latora V., Rapisarda A. and Ruffo 
S., Phys. Rev. Lett, 80 (1998) 692; Latora V., Rapis- 
arda A. and Tsallis C., Phys. Rev. E, 64 (2001) 056134; 
Pluchino A., Latora V. and Rapisarda A., Physica D, 
193 (2004) 315; Rapisarda A. and Pluchino A., Euro- 
physics News, 36 (2005) 202; Giansanti A., Moroni D. 
and Campa A., Physica A, 305 (2002) 137 ; J. Phys. A 
, 36 (2003) 6897 ; Chavanis P-H., Eur. J. Phys. B, 53 

(2006) 487; Antoniazzi A., Califano F., Fanelli D., 
and Ruffo S., Phys. Rev. Lett, 98 (2007) 150202; Morita 
H. and Kaneko K., Phys. Rev. Lett, 96 (2006) 050602; 

[3] Pluchino A. , Rapisarda A. and Tsallis C., Euro- 
physics Letters, 80 (2007) 26002 

[4] Pluchino A. , Rapisarda A. and Tsallis C., Physica 
A, 387 (2008) 3121 

[5] Tsallis C., J. Stat. Phys., 52 (1988) 479; Tsallis C., 
Gell-Mann M. and Sato Y., Europhys. News, 36 (2006) 
186, and references therein; Tsallis C., Milan J. Mathe- 
matics, 73 (2005) 145, and references therein 

[6] S. Umarov, C. Tsallis and S. Steinberg, cond- 
mat/0603593, Milan J. Math., 76 (2008) DOI 
10. 1007/s00032-008-0087-y 

[7] Tirnakli U., Beck C. and Tsallis C., Phys. Rev. E, 75 

(2007) 040106 (R); Tirnakli U., Tsallis C., Beck C., 
0802.1138, cond-matstat-mech, (2008) 

[8] F. Baldovin F. and A. Stella, Phys Rev. E, 75 (2007) 
020101 (R); Moyano L.G., Tsallis C., Gell-Mann M., 
Europhys. Lett, 73 (2006) 813; Thistleton W.J., 
Marsh J. A., Nelson K. and Tsallis C., IEEE Transac- 
tions on Information Theory , 53 (2007) 4805; Hilhorst 
H.J. and Schehr G., JSTAT, 06 (2007) P06003; ViG- 
nat C. and Plastino A., J. Phys. A , 40 (2007) F969; 
Rodriguez A., Schwammle V., Tsallis C., O8O4.1488 
[cond-mat.stat-mech] JSTAT, 2008 ( in press) 



